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Abstract The gravitational interaction between a protoplan- 
etary disc and planetary sized bodies that form within it 
leads to the exchange of angular momentum, resulting in 
migration of the planets and possible gap formation in the 
disc for more massive planets. In this article, we review the 
basic theory of disc-planet interactions, and discuss the re- 
sults of recent numerical simulations of planets embedded 
in protoplanetary discs. We consider the migration of low 
mass planets and recent developments in our understanding 
of so-called type I migration when a fuller treatment of the 
disc thermodynamics is included. We discuss the runaway 
migration of intermediate mass planets (so-called type III 
migration), and the migration of giant planets (type II mi- 
gration) and the associated gap formation in this disc. The 
availability of high performance computing facilities has en- 
abled global simulations of magnetised, turbulent discs to 
be computed, and we discuss recent results for both low and 
high mass planets embedded in such discs. 

Keywords planet formation • accretion discs • numerical 
simulations 



1 Introduction 

Planets are believed to form in the circumstellar discs of gas 
and dust that are observed around young stars (1). Within 
these discs, sub-micron sized interstellar dust grains col- 
lide and stick, slowly growing to form bodies comparable 
in mass to the Earth (Mq) via a multi-stage process which 
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involves bodies of increasing size colliding and sticking. 
These planetary bodies which form with masses up to ~ 10 
Me can accrete gas envelopes, leading to the formation of 
a gas giant planet (2), a process which, models predict, re- 
quires a few million years. The precursor ~ 10 Mq bodies 
which accrete gas are usual referred to as planetary cores, 
and their formation is most efficient at locations in the disc 
where the density of solids is highest. Thus gas giant plan- 
ets are thought to form in the cool disc regions beyond the 
'snowline', where volatiles can condense into ice grains. In 
most disc models, and apparently in our own early Solar 
System, the snowline is located at a radius ~ 3 AU from 
the central star The discovery of the first extrasolar planet 
in 1995 (3), a gas giant planet orbiting at 0.05 Astronomical 
Units (AU) from the star 5 1 Pegasi therefore came as a big 
surprise, and posed a problem for planet formation theory. 

The solution to this problem is thought to be that planets 
such as 5 1 Peg b probably did form several AU away from 
the star, but migrated inward afterwards, while still embed- 
ded in the circumstellar disc. The large number of extraso- 
lar planets now known (more than 350 at the time of writ- 
ing), with broad variation in their masses and orbital config- 
urations, indicates that migration is a common phenomenon 
during planet formation. 

Prior to the discovery of extrasolar planets, it was al- 
ready known that embedded bodies experience a force from 
the surrounding gas disc mm 

leading to orbital migra- 
tion. This can be appreciated in a qualitative way from Fig. 
[H where we show the surface density perturbation in a disc 
by a low-mass planet. The planet launches two tidal waves 
that are sheared into spiral waves by differential rotation in 
the disc. The gravitational pull from these waves leads to 
orbital migration of the planet. 

Disc-planet interaction has become the focus of intense 
research since 1995, and three types of migration can be dis- 
tinguished; Type I migration for low-mass planets (compara- 
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ble to the Earth), which do not strongly perturb the surround- 
ing disc. Type II migration for high-mass planets (compara- 
ble to Jupiter), which open deep gaps around their orbit, and 
Type III migration for intermediate-mass planets (compara- 
ble to Saturn), embedded in massive discs /s'). 

A significant fraction of the research in this area has 
utilised numerical simulations, and a number of insights have 
been gained from the results of such simulations which have 
subsequently been intepreted within the fraemwork of an an- 
alytic model. The increasing availability of parallel comput- 
ing facilities, combined with the development of simulation 
codes, has allowed increasingly complex simulations to be 
performed, including physics such as MHD turbulence and 
radiation transport. In this article we will review the basic 
theory of disc-planet interactions, highlight some recent de- 
velopments, and summarise the state of the field as it cur- 
rently stands. 

2 Type I migration 

Type I migration is thought to apply to low-mass planets 
(approximately 1-20 Mq). It is important to understand this 
regime, because not only does it apply to terrestrial planets, 
also the solid cores of gas giant planets have been subject 
to Type I migration before they obtained their gaseous enve- 
lope. This stage in giant planet formation requires more than 
10^ years (2). In this section, we will review the basic the- 
ory of Type I migration, focusing on laminar (non-turbulent) 
discs. 

Type I migration has long been regarded as the simplest 
migration mechanism to understand. Since low-mass plan- 
ets induce only small perturbations in the disc, a linear anal- 

§'s should provide accurate estimates of the migration rate 
. These can later be verified by numerical hydrodynami- 
cal simulations (e.g. 

2.1 Linear analysis 

It is convenient to express the perturbing potential due to the 
planet acting on the disc 0p as a Fourier series 

$>p(r,(p,f) = - = 'I>mir)cos{m(p-mQpt), (1) 

where (p is the azimuthal angle and is the orbital fre- 
quency of the planet, and to evaluate the response of the 
disc due to each component separately. For simplicity, we 
will consider only two-dimensional, isothermal discs in this 
section. The vertically integrated pressure is then given by 
P — c^E, where c is the sound speed and E is the surface 
density of the disc defined hy E ^ /"^ p dz- In vertical hy- 
drostatic equilibrium, the pressure scale height of the disc 
equals H — c/Q, where £2 is the angular velocity of the 



disc. A typical disc has h = H/r « 0.05. The unperturbed 
surface density follows a power law in radius, E °^ r^". 

When working with discs, angular momentum is a key 
quantity, dominating the dynamics. The angular momentum 
of the planet is Jp = MprpQp °^ y^r^ for a circular orbit. 
Therefore, changes in angular momentum are directly re- 
lated to a change in orbital radius, which makes the torque, 
the change in angular momentum, an important quantity to 
measure. 

Writing the perturbed surface density as 

E'{r,(l),t) ^E'„^{r)exp{im(p-imQpt), (2) 

one can numerically solve the linearised fluid equations to 
obtain the disc response. The torque on the planet, 

r = [ ErxVtppdS, (3) 

JOisc 

is directly related to the imaginary part of the surface density 
perturbation: 

^ = -nm0,„J^m{E'J, (4) 

and the torque r,„ can be found by integrating over the whole 
disc. The total torque is equal to the sum over all I^,. 

Neglecting pressure effects, the disc can exchange an- 
gular momentum with the planet at a Lindblad resonance, 
which is located at a position in the disc where m{£2 —£2p) = 
zLk, or at a corotation resonance, where Q — Qp. Here, K is 
the epicyclic frequency, equal to i2 in a Keplerian disc. At 
a Lindblad resonance, a pressure wave is excited that propa- 
gates away from the planet, which cause an oscillatory disc 
response in the solutions to linearised equations described 
above. The waves transport angular momentum, resulting in 
a torque on the planet. The outer wave removes angular mo- 
mentum from the planet, while the inner wave gives angular 
momentum to the planet. The torque therefore depends on 
the difference in the strengths of the two waves. 

There are several effects that conspire to make the outer 
wave stronger , the most important one being that outer 
Lindblad resonances are located closer to the planet than in- 
ner Lindblad resonances. Therefore, for all reasonable sur- 
face density profiles the total Lindblad torque on the planet 
is negative, resulting in inward migration. The most recent 
linear calculations for a 2D disc (9) resulted in a Lindblad 
torque 

rL.ii„ = - (3.2 + 1.468a) |2^p'-p^2p', (5) 

where q denotes the mass of the planet in units of the central 
mass, rp denotes the orbital radius of the planet and Ep is the 
surface density at the location of the planet. Note that a is 
usually thought to be positive, resulting in a negative torque 
on the planet. 
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Fig. 1 Surface density perturbation for a 4 planet (located at 
{x,y) = (— 1,0)) embedded in a disc with h = 0.05, showing the promi- 
nent spiral wakes associated with Lindblad torques. 



In figure[T]we show the surface density perturbation due 
to a low-mass planet from a numerical hydrodynamical cal- 
culation. Clearly visible are the two spiral waves launched 
by the planet that result in the above Lindblad torque. 

As mentioned above, angular momentum exchange also 
takes place at a corotation resonance. Since the location of 
this resonance is very close to the orbit of the planet, the 
resulting corotation torque is potentially very strong. Since 
the relative velocity of the gas and the planet is subsonic 
near the corotation resonance, the planet is unable to excite 
waves at this location ClZ) . Semi-analytical calculations 1^ 
result in a linear corotation torque 



'c.lin 



1.36 



and therefore a total linear torque 



L.lin 



c.lin 



(1.16 + 2.828a)^i;pr^i2^. 



(6) 



(7) 



Again, since a is usually thought to be positive, we find that 
low-mass planets will migrate inward. In equation |6] the 
quantity in brackets is the power law index of the specific 
vorticity, or vortensity (V x \)/E, where v is the velocity 
in the disc. The corresponding migration rate is (assuming a 
circular orbit): 



=- (1.16 + 2.828a) 4 



rpX2pMp 



p^'p. 



(8) 



A major problem with the resulting inward migration 
rates is that they are much too high: Type I migration would 
quickly take all low-mass planets, including the cores of gas 
giant planets, very close to or even into the central star (11 ih . 



Models of planet formation including Type I migration (Il3l) 
have great difficulty explaining the observed distribution of 
planets. 

Several mechanisms have been proposed to slow down 
Type I migration. Within linear theory, it is first of all im- 
portant to note that the two-dimensional approximation is 
not valid once the gravitational sphere of influence of the 
planet (the Hill sphere) is smaller than the pressure scale 
height of the disc. It is easy to see that for a small enough 
planet, the upper layers of the disc do not feel the gravita- 
tional pull of the planet, while in the above two-dimensional 
theory they fully contribute to the torque. A torque formula 
that accounts for this effect reads (9) 

n,2 



r^DM = - (1.364 + 0.541a) |2 2:pr;i2p^ 

Note that the difference between the 2D and the 3D result 
strongly depends on the surface density profile. In favourable 
cases, migration speeds can be reduced by a factor of a few. 
In the 2D numerical simulations presented below, a gravita- 
tional softening parameter b is used to account for 3D effects 
in an approximate way: 



(9) 



GMr, 



(10) 



It is expected that for b ^ h, 3D effects on the Lindblad 
torque can be captured approximately. The same may not be 
true for the non-linear torque discussed below. In (14), equa- 
tion|9]has been successfully compared to 3D isothermal nu- 
merical simulations. Note, however, that their adopted den- 
sity profile allows for weak corotation torques only. 

The above analysis was simplified by neglecting mag- 
netic fields, self-gravity and detailed thermodynamics. In- 
cluding magnetic resonances in the linear analysis can re- 
duce the torque significantly if a toroidal magnetic field is 
present (15). On the other hand, self-gravity tends to accel- 
erate Type I migration due to a shift in resonances (3 17). 
More recently, releasing the assumption of an isothermal 
disc has been the subject of intensive study, since simula- 
tions presented by (18) which included radiation transport 
in the disc suggest that type I migration may be stopped or 
reversed in a disc which cools inefficiently. A linear analysis 
([19') suggested that a linear effect due to the presence of a ra- 
dial entropy gradient may reduce or even reverse migration. 
However, it was subsequently shown (|20|) that this effect is 
in fact small and can not account for the results obtained in 



2.2 Beyond the linear model 

Although intuitively, one would expect that low-mass plan- 
ets induce linear perturbations only, recent work has shown 
that this is not the case. In a sense, introducing a planet in 
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Fig. 2 Schematic view of the streamline pattern close to the planet, 
which is located at {r,(p) = (1, 7t) and is denoted by the black circle. 
Material inside the black streamline (the separatrix) executes horseshoe 
turns on both sides of the planet. 
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Fig. 3 Torques on a = 1 .26 • 10^' planet (4 around a Solar mass 
star) embedded in an inviscid disc with h = 0.05. All torques have been 
normalised to Zp = rp = X2p = 1 . Results for different density profiles 
are shown, with dotted horizontal lines indicating results obtained from 
linear theory. A gravitational softening of b/h = 0.6 was used. 



a disc amounts to a singular perturbation at corotation, sug- 
gesting that the corotation torque can differ from its hnear 
value. 

The corotation torque is due to material that, on average, 
corotates with the planet. Two ways of obtaining the corota- 
tion torque can be found in the literature: one can perform a 
linear analysis of perturbed circular orbits around the coro- 
tation resonance (see above), or, alternatively, one can 
look at the expected pattern of the streamlines close to the 
planet (see figure m and analyse the torque due to material 
executing horseshoe turns (12 lb . We will denote the former 
result the linear corotation torque, and the latter the horse- 
shoe drag. It is important to stress that there is no horseshoe 
region in linear theory, since these bends are a genuinely 
non-linear phenomenon. 

2.2.1 Isothermal results 

Considering again two-dimensional, isothermal (or, more gen- 
eral, barotropic) flow, in which case specific vorticity is con- 
served along streamlines, it is easy to see that when execut- 
ing the turn, material should change its density. Since it is 
forced to a different orbit by the planet, the vorticity of the 
material changes, and conservation of vortensity then dic- 
tates that there should be an associated change in density. 
This change in density is of opposite sign on the different 
sides of the planet, hence the planet feels a torque. This 
torque depends sensitively on the width of the horseshoe re- 
gion .Xs U^- 

Note that Xs is in units of r^ . Just as in the linear case, the 
horseshoe drag is proportional to the radial gradient in spe- 
cific vorticity. In the special case of a = 3/2, the density 



change needed to compensate for the change in vorticity 
exactly amounts to the background profile, resulting in no 
torque on the planet. 

Until recently, it has never been clarified how equations 
|6]and[TT]are related. There is no reason why the value of Xs, 
unspecified in the original theory (|2ll). should adjust itself 
to match the linear torque. Detailed modeling of the horse- 
shoe region shows that Xs °= \fqjh, so that the horse- 
shoe drag scales in exactly the same way as the linear coro- 
tation torque. The only way they differ is in magnitude, and 
in general the horseshoe drag is much stronger (23). 

It is intuitively clear that these two distinct torques can 
not exist at the same time: material either follows a per- 
turbed circular orbit or executes a horseshoe turn. At suf- 
ficiently early times after inserting the planet in the disc, 
therefore, before any turns can have occurred, one would 
expect linear theory to be valid. Once material starts to exe- 
cute horseshoe turns, the linear corotation torque will get re- 
placed by the non-linear horseshoe drag, which, as remarked 
above, will in general be much stronger. 

This can be verified using numerical hydrodynamic sim- 
ulations. Resolution is crucial, since for low-mass planets 
Xs <C /i <C 1, and in order to obtain reliable estimates of the 
torque, the horseshoe bends need to be resolved. Simulation 
results for different values of a are displayed in figure [3] 
showing the total torque on the planet as a function of time. 
Overplotted are the results from solving the linear equations 
(dotted lines). It is clear that only in the case where o: = 3/2 
the non-linear hydrodynamic calculations match linear the- 
ory. In this special case, the corotation torque vanishes. For 
all other values of a, a departure from linear theory can be 
seen, pointing to a stronger corotation torque than expected 
from linear theory. However, the total torque is consistent 
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with the linear corotation torque replaced by the stronger 
horseshoe drag, showing that indeed non-linear effects play 
an important role in the corotation region. 

Another important aspect is the time scale to set up the 
torque. Since the linear corotation torque has no definite 
scale, it takes approximately one orbit to develop. The horse- 
shoe drag, on the other hand, does have an intrinsic scale, 
namely Xj. An associated time scale is the libration time, 
which is basically the time it takes for a fluid element, ini- 
tially at a distance Xs from the orbit of the planet, to complete 
one whole horseshoe orbit: 

Tiib = :r-r2p-i. (12) 

It takes a fraction of the libration time for the horseshoe drag 
to develop, which is therefore different for different planet 
masses since Xs °^ ^Jq. The transition between the linear re- 
sult (obtained after approximately 1-2 orbits), and the more 
slowly developing non-linear horseshoe drag can be clearly 
seen in figure [3] 
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Fig. 4 Torques on a = 1 .26 • 10^^ planet (4 M^i around a Solar mass 
star) embedded in an inviscid disc with h = 0.05. All torques have 
been normalised to Zp = rp = i2p = 1. The black line indicates the 
isothermal result, grey curves results from adiabatic calculations for 
two different temperature profiles. A gravitational softening of bjh = 
0.4 was used. 



2.2.2 Extension to non-isothennal discs 

These simple models of Type 1 migration lack a crucial in- 
gredient: a proper treatment of the disc thermodynamics. 
Most previous simulations were indeed done in the locally 
isothermal limit, where the temperature of the disc is a pre- 
scribed and fixed function of radius. Such a model would 
apply if any excess energy could be radiated away very ef- 
ficiently. However, in discs that can not cool efficiently (i.e. 
are optically thick), the presence of a radial entropy gradient 
leads to strong coorbital torques that can reverse the direc- 
tion of migration [l^. This can again be understood in terms 
of the horseshoe drag. When the gas is not able to cool (the 
adiabatic limit), entropy is conserved along streamlines. In 
a similar way as the vortensity gradient led to the barotropic 
horseshoe drag, a radial entropy gradient leads to a strong 
horseshoe drag related to density changes as the fluid tries to 
conserve its entropy when making the horseshoe turn (j20t). 
When the entropy decreases outward, this horseshoe drag 
gives a positive contribution to the torque. 

This is illustrated in figure |4] where we compare adia- 
batic and isothermal simulations with a = 3/2 for different 
temperature profiles (parametrised by j3 = —dlogT/dlogr). 
An adiabatic exponent 7=5/3 was chosen. Due to the isen- 
tropic nature of linear waves, the density of the wakes is 
lowered by a factor y compared to isothermal simulations, 
which makes the Lindblad torque a factor y smaller. Note 
that the density profile is such that the barotropic corotation 
torque vanishes. The entropy follows a power law, initially, 
with index — ^ = —j5 + (7— l)a. Therefore, the isothermal 
disc (/3 = 0) with no vortensity gradient (a = 3/2) does have 
an entropy gradient (^ = — 1), giving rise to a strong horse- 
shoe drag. This is clear from the grey solid line, which shows 
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Fig. 5 Torques on a = 1 .26 • 10^^ planet (4 Mca around a Solar mass 
star) embedded in a 3D inviscid disc with /; = 0.05. All torques have 
been normalised to Zp = rp = i2p = 1. Results for different midplane 
densities are shown, and the differences between the curves are a re- 
sult of the corresponding change in opacity, which affects the cooling 
properties of the disc. 



a large swing between 2 and 10 orbits as a result. When we 
modify the temperature structure to reverse the entropy gra- 
dient, a corresponding opposite swing is observed (dashed 
grey line in figure |4]i. Even without an additional contribu- 
tion from a vortensity gradient, the torque is reduced by a 
factor of 3 compared to the linear value. 

Finally, we release the two-dimensional approximation, 
and enable the disc to cool through radiation. The resulting 
3D radiation-hydrodynamical models (llSh are the most re- 
alistic simulations of low-mass planet migration so far In 
figure |5] we show the results for a = 1/2 and )3 = 1. Dif- 
ferent midplane densities Pp were used, where Pp — 10^^' 
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g cm^^ approximately corresponds to the minimum mass 
Solar nebula at 5 AU. Note that the torques are still nor- 
malised to i^p = 1 . Therefore, different densities mainly in- 
dicate different opacities. 

The high density case is almost equivalent to an adia- 
batic simulation. This is easy to understand, since the disc 
is very optically thick in this case and basically can not cool 
through radiation. This again leads to a positive torque, and 
outward migration. Lowering the density leads to more effi- 
cient cooling, and one would expect to return to the isother- 
mal result when the opacity is low enough. This is indeed 
the case: upon decreasing the density by a factor of 10 the 
torque is already negative, and with another factor of 10 
we can almost reproduce the isothermal result. Therefore, 
the entropy-related horseshoe drag is important in the dense 
inner regions of protoplanetary discs. This mechanism can 
then act as a safety net for low-mass planets, since they stop 
their inward journey as soon as they reach the dense inner 
regions of the disc. 



2.2.3 Saturation 



2.4 Concluding remarks 

Type I planetary migration is not as well-understood as pre- 
viously thought. Even in the simple isothermal case, strong 
deviations from linear theory occur whenever the disc has a 
radial gradient in specific vorticity. In other words, the coro- 
tation torque is always non-linear. Therefore, one should be 
careful in applying the widely used linear torque formula 
presented in (9)- Even more so for non-isothermal discs, 
where a radial entropy gradient can contribute to the non- 
linear torque to the effect of reversing the direction of migra- 
tion. In order for the torque to be sustained, some diffusion 
of energy and momentum is needed. 



3 Type III migration 

The most recently discovered migration mechanism, first re- 
ferred to as runaway migration (|8|), relies on strong corota- 
tion torques induced by the radial movement of the planet 
itself as it migrates. Later, it has become known as Type 111 
migration (e.g.uSh. 



In absence of any diffusive process in the disc, the horse- 
shoe region is a closed system, and it can therefore only give 
a finite amount of angular momentum to the planet. Phase 
mixing leads to a flattening of the vortensity/entropy pro- 
files, which results in a vanishing torque after the gas has 
made several horseshoe turns. In other words, the corota- 



tion torque saturates (1251) . The horseshoe region needs to be 
fed fresh material in order to sustain the torque. This can be 
done through viscous diffusion in the isothermal case ^2^}, 
and it has been observed that by including a small kinematic 
viscosity and heat diffusion the torque can be sustained (l2ol) . 
This has also been shown with radiative cooling instead of 
heat diffusion in illh . 



2.3 Outstanding issues 

The non-linear torque strongly depends on the width of the 
horseshoe region Xg. It is therefore of critical importance 
to have a good physical understanding of what determines 
this width when three-dimensional effects are taken into ac- 
count. 

It is not clear yet how the entropy and vortensity related 
torques work together. Vortensity is not conserved in general 
along streamlines in 3D non-barotropic discs, but neverthe- 
less there appears to be a contribution from the radial gradi- 
ent in vortensity to the horseshoe drag in this case. A proper 
understanding of these processes will lead to an improved 
migration law for low-mass planets. 



3.1 Theory 

We will return to the simple case of a laminar, isothermal 
disc. Consider the horseshoe region as depicted in the left 
panel of figure |6] and for now assume that it has constant 
specific vorticity, either through saturation or due to the ini- 
tial density profile. The corotation torque therefore vanishes. 
If we now let the planet migrate radially, the horseshoe re- 
gion as viewed in a frame comoving with the planet changes 
its shape (middle panel of figure |6]l. If the planet is moving 
inward, as is the case in figure |6] some fresh material from 
the inner disc will enter the horseshoe region, make the turn 
and move into the outer disc, exerting a torque on the planet 
in the process. It is easy to see that this torque is proportional 
to the migration speed, and acts to accelerate the migration. 
There is therefore a positive feedback loop, and the possibil- 
ity of a runaway process. 

The mass flow relative to a planet migrating at a rate rp 
(we assume a circular orbit, so that the semimajor axis a = 
rp) is given by InL^rp, with the surface density at the up- 
stream separatrix. Since we assume the ordinary, vortensity- 
related corotation torque to be fully saturated, the only horse- 
shoe drag on the planet comes from this fresh material exe- 
cuting a horseshoe turn, giving rise to a torque ( 1291) : 



r^,3 = InZsi-pQprpXs 



(13) 



The system of interest that is migrating is the planet itself, 
mass of gas bound within its Roche lobe, and the horseshoe 
region. The drift rate of this system is given by 



(Mp+MHs+MR)i2| 



.p'"p'"p 



4KrlxAQ„rr,rp + IF^, (14) 



7 




1 








'-_ 4 - / 








H 1 -/ 




: oil 





Fig. 6 Streamlines near the corotation region for different migration speeds, in a frame moving with the planet. Left panel: no migration, middle 
panel: slow inward migration, right panel: fast inward migration. 
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Fig. 7 Non-dimensional migration rate as a function of the mass deficit 
parameter . 



which can be rewritten as 
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Fig. 8 The three migration regimes in disc-planet mass parameter 
space, for h = 0.05. 



where denotes the coorbital mass deficit divided by the 
planet mass, 



Mo ' 



2ri_ 



(15) 



where we have redefined Mp as the mass of the planet plus 
the mass within its Roche lobe. Here, II is the Lindblad 
torque, which may deviate from the linear Lindblad torque 
due to shock formation of the appearance of a gap around 
the orbit of the planetl. The numerator of the first term on 
the right hand side is called the coorbital mass deficit 
Basically, this is the difference between the mass the horse- 
shoe region would have in a disc with Z = l^s and its ac- 
tual mass. It is clear that when this coorbital mass deficit 
becomes comparable to Mp, very high migration rates are 
possible. However, in this limit the above analysis becomes 
invalid. A more careful treatment (.29; i30i) results in: 



-M4sign(Z)(l-(l-|Z|<i)'/') = 



2rL 



Mpi2prprp^f 



(16) 



34 



(17) 



is the drift rate required for the planet to migrate a distance 
Xs in one libration time, and Z = rp/r^.f- A quantity with sub- 
script < 1 indicates that the minimum of that quantity and 1 
is to be taken. Equation[l6]gives rise to interesting behavior. 
Assuming Lindblad torques to be small, solutions to equa- 
tion[T6]are plotted in figure|2l ForM^ < 1, only one solution, 
Z = 0, is possible. There is a bifurcation at M^ = 1, beyond 
which there are three possible solutions: steady inward or 
outward Type III migration, plus an unstable solution Z = 0. 
Incorporating 11 in this picture removes this symmetry be- 
tween inward and outward migration {3^), favouring migra- 
tion towards the central star, although outward migration is 
still possible (.31.) . 
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Fig. 9 This figure shows a snapshot of a planet with inital mass of 1 Jupiter mass embedded in a protoplanetary disc. Note the formation of non 
linear spiral waves and the accompanying gap. 



3.2 Occurence 

There are different ways in which a coorbital mass deficit 
can be obtained. The first possibility is when the planet opens 
up a (partial) gap (8). Clearly, in this case there is a dif- 
ference between the actual mass of the horseshoe region 
and the mass of the horseshoe region with E — Eg, hence 
there is a coorbital mass deficit. For a high enough disc 
mass, this deficit can become comparable to the planet mass. 
More massive planets require a more massive disc to en- 
ter the Type III migration regime. On the other hand, low- 
mass planets do not open up a gap and therefore their coor- 
bital mass deficit is zero. It turns out that in a typical disc, 
Saturn-mass planets are most susceptible to Type III mi- 
gration (8). This is depicted in figure [8] where we show 
the regime of Type III migration for different disc masses. 
Here, /id = nr^Ep/M^, is a dimensionless disc mass param- 
eter, and Q is the corresponding Toomre stability parameter 
against fragmentation. For the Minimum Mass Solar Nebula 
(MMSN) and h = 0.05, g > 10 inside 10 AU, making Type 
III migration unlikely. However, the MMSN is effectively 
a lower limit on the true disc mass profile, and for massive 
discs Type III migration may be an important effect. 

Another way for a coorbital mass deficit to arise is when 
the background disc has a sharp edge. This was studied for 
Jupiter-mass planets (q = 0.001) in ((28; 31). Such a coor- 
bital mass deficit would also affect low-mass planets, but 
this has not been studied so far. 

Type III migration is a numerically challenging problem, 
since it strongly depends on material close to the planet. For 
example, in the above analysis we redefined the planet mass 
to include all material inside the Roche lobe. If the planet ob- 
tains a significant envelope, torques from within the Roche 



lobe can be large (32) and oppose Type III migration. How- 
ever, since self-gravity is neglected, this material has effec- 
tively no inertia, making this calculation inconsistent ( 1291) . 

To make further progress and deal with such a massive 
envelope, it is necessary to incorporate effects of self-gravity, 
at least in an approximate way (33). Moreover, to avoid an 
artificially large mass build-up inside the Roche lobe due to 
the isothermal approximation, a temperature profile for the 
envelope needs to be supplied. Only then can simulations of 
Type III migration, including torques from inside the Roche 
lobe, be shown to be converged (.33). 

Type III migration can be very fast. The typical time 
scale is '"p/'"p,f = 2/3^ orbits. For — h, appropriate for 
Saturn-mass planets 1221), the migration timescale is approx- 
imately 200 orbits, or 2360 year for a planet starting at 5.2 
AU. However, it is not clear at this point how long the coor- 
bital mass deficit can be retained. For example, unless the 
surface density profile decreases steeply with radius, the coor- 
bital mass deficit will naturally decrease as the planet moves 
inward because of the changing factor r^E^ in equation [Tsl 
Outward Type III migration does not suffer from this effect, 
but simulations show that in this case the planet acquires a 
very massive envelope that pulls it into the Type II regime 
^1). Therefore, it seems that Type III migration is of limited 
importance for high-mass planets, changing the semimajor 
axis by a factor of a few (l28tl3lb . 



4 Type II migration in viscous, laminar discs 

When a planet grows in mass, the spiral density waves ex- 
cited by the planet in the disc can no longer be treated as 
linear perturbations in the manner described in section |2l 
The planetary wake becomes a shock-wave in the vicinity 
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Fig. 10 The left panel shows the orbital radius of giant planet as it migrates in the disc model described in the text. The right panel shows the 
mass of the planet as it concurrently accretes gas from the disc. 



of the Lindblad resonances where it is launched. Shock dis- 
sipation, as well as the action of viscosity, leads to the depo- 
sition of angular momentum associated with the spiral wave 
locally in the disc. Material is thus pushed away from the lo- 
cation of the planet on either side of its orbital location, and 
an annular gap begins to form, centred on the planet semi- 
major axis. The equilibrium structure and width of the gap 
is determined by a balance between the gap-closing viscous 
and pressure forces and gap-opening gravitational torques 
(e.g. 



34) 



In order for a gap to form and be maintained, two ba- 
sic conditions must be satisfied. The first is that the disc 
response, via the launching of spiral wakes, should be non 
linear, such that gravitational forces due to the planet over- 
whelm pressure forces in the neighbourhood of the planet. 
This condition for non linearity is equivalent to the planet 
Hill sphere radius exceeding the disc vertical scale height, 
and may be expressed mathematically as: 



Mp 

3Mh. 



1/3 



>H. 



(18) 



where a is the planet semimajor axis. The second condition 
is that angular momentum transport by viscous stresses in 
the disc be smaller than planetary tidal forces. This can be 
expressed approximately as 



> 



(19) 



V 

where v is the kinematic viscosity. For more discussion about 
these aspects of gap opening, see the following papers (35; 

SHIii). 

For standard disc parameters, H/r — 0.05, and v / {a^Q ) = 
10"^, equation ([18]) givesMp/M* > 3.75 x 10""*, while equa- 



with a mass exceeding that of Saturn will begin to open a 
noticeable gap. 

Because the disc response to a massive embedded pro- 
toplanet is non linear, the study of the orbital evolution and 
gas accretion by such a planet has been the subject of exten- 
sive study via numerical simulations. Both 2D (e.g^36j; 38; 
37; 39) and 3D simulations jiol: Wlh have been performed, 
and overall the results found in these studies are in good 
agreement. Differences obviously arise in the details of the 
gas flow near the planet when comparing 2D and 3D simu- 
lations, but the global properties of the models are otherwise 
similar as the gas is pushed away from the planet vicinity by 
the process of gap opening. 

Two snapshots from a 2D simulation of a disc with H/r = 
0.05 and v/{a^Q) — 10"^ ai'e shown in figure |9] where 
the initial mass of the planet Mp — 1 Mjup (where Mjup is 
Jupiter's mass). The formation of a deep annular gap in the 
vicinity of the planet is apparent. The parameters used in this 
simulation are such that both the condition for non linearity 
and the condition for tides to overwhelm viscous stresses are 
satisfied, as described above. 

The evolution of the planet semimajor axis is shown in 
the left panel of figure [TO] As the gap forms and the planet 
migrates inward, the migration rate drops significantly be- 
low that which is predicted for type I migration. Once a 
clean gap is formed, the planet orbital evolution is controlled 
by the rate at which the outer gap edge viscously evolves, 
such that the migration time, Tmig, can be approximated by 



3v ' 



(20) 



For typical disc parameters the time to migrate from a dis- 



tion ( fT9] ) gives Mp/M* > 4 x 10 The estimates are in rea- tance of rp = 5 AU to the central star is ~ 2 x 10^ years. 



sonable agreement with the results of numerical simulations 
of gap opening by (36) and ( i38i) . Accordingly, for typical 
protoplanetary disk parameters, we can expect that a planet 



As the planet mass increases, or the disc mass decreases, 
the inertia of the planet can become important and begin to 
slow the migration. This occurs when the local disc mass 
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contained within the orbital radius of the planet falls beneath 
the planet mass (7Zr^E{rp) < Mp). In this case the migra- 
tion rate is controlled by the rate at which mass can build- 
up at the edge of the gap such that the disc mass becomes 
comparable to the planet mass locally The slowing 
of migration is apparent in the left panel of figure [TOl and 
a simple linear extrapolation of the migration history indi- 
cates that the time required for this planet to migrate into 
the central star is ~ 3 x 10"* orbits (or ~ 3 x 10^ years for 
a planet initially located at 5 AU. The right hand panel of 
figure [To] shows the evolution of the planet mass as it ac- 
cretes gas from the disc. The mass accretion rate is slowed 
by the formation of the gap, but nonetheless gas is still able 
to diffuse through the gap onto the planet allowing it to 
grow. Linear extrapolation of the planet mass shown in fig- 
ure [To] indicates that the planet will have grown to a mass 
of Mp ~ 3.5Mjup by the time it reaches the central star, such 
that migration and concurrent gas accretion occur during the 
gap opening phase of giant planet formation. 

The simulations for gap forming planets presented in 
this section indicate that planets grow substantially through 
continued gas accretion as they migrate. This raises an inter- 
esting question about the type of correlation which expected 
between planetary semimajor axes and masses. Taken at face 
value, a model which predicts that planets grow as they mi- 
grate inward will predict an inverse correlation between plan- 
etary mass and semimajor axes. The strength of this cor- 
relation will be weakened by variations in protoplanetary 
disc masses, planet formation times, disc lifetimes etc, but 
if these parameters are peaked around certain typical values 
in nature, then the correlation should be apparent in the ex- 
trasolar planet data. The absence of such a correlation may 
in part be explained by post-formation interaction within a 
planetary system, leading to significant gravitational scat- 
tering. The possibility that significant migration may occur 
prior to the last stages of giant planet formation may also 
be important. Nonetheless, indirect evidence for the past oc- 
curence of Type II migration in some extrasolar planetary 
systems is provided by multi-planet systems which display 
the 2:1 mean motion resonance (or other resonances of low 
degree). Such a configuration is most naturally explained by 
two giant planets undergoing convergent migration at rates 
predicted by Type II migration theory, as illustrated by the 
numerical simulations presented in (e.g.l43l: l44l). 



important effects (l45l) . More work is necessary in this area, 
also to completely release the isothermal assumption to ob- 
tain a reahstic envelope structure for the planet. 



5 Turbulent protoplanetary discs 

The majority of calculations examining the interaction be- 
tween protoplanets and protoplanetary discs have assumed 
that the disc is laminar The mass accretion rates inferred 
from observations of young stars, however, require an anoma- 
lous source of viscosity to operate in these discs. The most 
likely source of angular momentum transport is magnetohy- 
drodynamic turbulence generated by the magnetorotational 
instability (MRI) (46) Numerical simulations have been per- 
formed using the local shearing box approximation (see ref- 
erence '4?', for a review) and using global geometry (e.g. [4^. 
In each of these set-ups, runs have been conducted using non 
stratified models in which the vertical component of grav- 
ity is ignored, and in more recent work the vertical compo- 
nent of gravity has been included in shearing box (.49i) and 
global models (l50i and references therein). These model in- 
dicate that the nonlinear outcome of the MRI is vigorous 
turbulence, and dynamo action. Recent work has indicated 
that the picture is more complicated than previously thought, 
with there being a strong dependence of the saturated state 
of the turbulence on the value of the magnetic Prandtl num- 
ber(51;^. 




4. 1 Outstanding issues 

It has turned out to be difficult to improve on the non-self- 
gravitating, isothermal results in which gas inside the Roche 
lobe is considered to be part of the planet ( 8). Approximate 
corrections for self-gravity and non-isothermality are neces- 
sary to make progress on high-mass planets (33). The first 
fully self-gravitating calculations indicate that this can have 



Fig. 11 This figure shows a snapshot of a 10 Ma, protoplanets embed- 
ded in the turbulent protoplanetary disc model described in the text. 

Prior to the reaUsation that the strength of the MHD tur- 
bulence depends on the Prandtl number, it was recongnised 
that non ideal MHD effects might be important in proto- 
planetary discs. The ionisation fraction in cool, dense pro- 
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toplanetary discs is expected to be small in the planet form- 
ing region between 0.5 - 10 AU. Only the surface layers 
of the disc, which are exposed to external sources of ioni- 
sation such as X-rays from the central star or cosmic rays, 
are likely to be suffici ently ionised to sustain MHD turbu- 
lence (e.g. i53i i54j; |55|; l56h . Predictions about the structure 
and depth of the interior non turbulent "dead zone" depend 
on complex chemical reaction networks and the degree of 
dust grain depletion, issues which themselves may be in- 
fluenced by the process of planet formation via dust grain 
growth. 

Although there are a number of uncertaintities remaining 
about the nature and strength of disc turbulence, it is clear 
that an accurate understanding of planet formation and disc- 
planet interaction requires the inclusion of turbulence in the 
models. Work conducted over the last four or five years has 
examined the interaction between planets of various masses 
and turbulent protoplanetary discs, where these studies have 
usually simulated explicitly MHD turbulence arising from 
the MRI. We now review the results of these studies, and 
present some more recent models. 



Semimajor axis versus time 
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Fig. 13 This shows the variation of semimajor axis with time (mea- 
sured in planet orbits at r = 3). The green line shows the planet evolu- 
tion in a turbulent disc, the red line shows the evolution in an equivalent 
laminar disc. 



5.1 Low mass protoplanets in turbulent discs 

The interaction between low mass, non gap forming proto- 
glanets and turbulent discs has been examined by ( 10,; JT; 
58t 1591) These calculations show that interaction between 
embedded planets and density fluctuations generated by MHD 
turbulence can significantly modify type I migration, at least 
over time scales equal to the duration of simulations that are 
currently feasible (f ^ 150 planet orbits). The influence of 
the turbulence is to induce a process of 'stochastic migra- 
tion' rather than the monotonic inward drift expected for 



planets in laminar discs. The simulations presented in the 
papers cited above considered the interaction of planets with 
non stratified (no vertical component of gravity) disc mod- 
els (for reasons of computational expense), but the work of 
indicates that the amplitude of midplane density fluctu- 
ations in stratified models is weaker than in their non strat- 
ified counterparts because of magnetic buoyancy effects re- 
moving net magnetic flux from the simulation domain. It is 
therefore important to consider the evolution of planets in 
stratified models. 

One such model is presented in figure [TT] which shows 
a snapshot of the midplane density in a turbulent disc with 
an embedded 10 Earth mass planet. The disc model in this 
case has aspect ratio H/r = 0.1, and was initiated with a ra- 
tio of magnetic to thermal gas pressure Pg/Pm = 25 though- 
out the computational domain. As such the model is very 
similar to those presented in ( fsol) . The disc extends from 
an inner radius of 1.6 AU to 16.6 AU, and has a vertical 
domain covering a total of 9 scale heights. The resolution 
corresponds to (A^^, A^^., A^e) = (464, 280, 1200). Figure [m 
shows the turbulent density fluctuations are of higher am- 
plitude than the spiral wakes generated by the planet, as the 
planet induced wakes are not visible against the perturba- 
tions induced by the MHD turbulence (see alsoflOtlssI). In- 
deed, density fluctuations generated by turbulence in simu- 
lations are typically 5E/L ~ 0.15 - 0.3, with peak fluctua- 
tions being ^(1). Thus, on a length scale equivalent to the 
disc thickness H, density fluctuations can contain more than 
an Earth mass in a disc model a few times more massive than 
a minimum mass nebula. 

The left panel of figure [12] shows the variation in the 
torque per unit mass experienced by the planet shown in fig- 
ure [TT] The green line shows the torque experienced by a 
planet in the turbulent disc model, and the red line shows 
the torque experienced by the same planet in an equivalent 
laminar disc model (the torque value here is ~ 3 x 10^^). 
It is clear that the planet in the turbulent disc experiences 
strongly varying stochastic torques whose magnitudes are 
significantly larger than the expected type I torques. The 
right hand panel shows the running time average of the torque 
experienced by the planet. The light blue, green and dark 
blue lines correspond to the torque from the outer disc, the 
inner disc, and the total torque, respectively. The red lines 
show the equivalent qualities from the laminar disc model. 
It is clear that while the various torque contributions are of 
similar magnitudes when comparing the laminar and turbu- 
lent models, the net average torque in the turbulent model 
is positive, indicating that the planet has migrated outward 
slightly over the duration of this simulation. 

This expectation is confirmed by figure [131 which shows 
planet semimajor axes for both turbulent (green line) and 
laminar (red line) models. It is clear that for the duration 
of this run the planet in the turbulent disc has undergone 
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Fig. 12 The green line in left panel shows the torque experienced hy the planet emhedded in the turbulent disc. The red line shows the torque 
experienced by a planet embedded in an equivalent laminar disc. The right panel shows the time averaged values of the torques. From top to 
bottom, red lines are the torques from the laminar disc calculation corresponding to the disc interior to the planet, the total torque, and the disc 
exterior to the planet. Similarly the turquoise, purpole and green lines are from the turbulent disc simulation. 
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Fig. 14 The left panel shows the distribution of torques experienced by the planet embedded in the turbulent disc. The right panel shows the 
autocorrelation function of the torque plotted in the left panel of Fig. ll3l 



Stochastic migration, resulting in it having moved outward 
sUghtly. A key question is whether the stochastic torques can 
continue to overcome type I migration over time scales up 
to disc life times. A definitive answer will require very long 
global simulations, but simple estimates can be made based 
on the simulated torques and their statistical properties. 

If we assume that the planet experiences type I torques 
from the turbulent disc, on top of which are linearly super- 
posed Gaussian distributed fluctuations with a characteristic 
correlation time Tc, then we can write an expression for the 
time averaged torque experienced by the protoplanet, T, in 
the form 

T =<T > (21) 

v'tot 

where Oj is the standard deviation of the torque amplitude, 
< r > is the underlying type I torque, and ftot is the to- 
tal time elapsed, measured in units of the torque correlation 
time Tc. In other words we treat the calculation as a simple 
signal-to-noise problem. Convergence toward the underly- 
ing type I value is expected to begin once the two terms on 



the right hand side become equal, as the fluctuating com- 
ponent should begin to average out. The torque distribution 
function for this simulation is shown in the left panel of fig- 
ure [141 and the torque autocorrelation function is shown in 
the right hand panel. Expressed in code units, the rms torque 
value is ^ 4 X 10^^ and the correlation time is Xmathrmc ~ 1 
planet orbital period. The type I torque is < T >~ 3.5 x 
10"^, so the time over which type I torques should be appar- 
ent is 130 orbits, slightly more than the simulation run time. 
This prediction suggests that we should expect to see the 
planet begin to migrate inward in figure [13] but this clearly 
is not the case. 

The lack of inward migration in figure [13] may just be 
a consequence of not having run the simulation for long 
enough, such that the fluctuating torques have not had time 
to begin to average out. The similarity between the run time 
of the simulation and the estimate of the time when type I 
effects should become apparent suggests that this may be 
the case. If we look at the right panel of figure [12] how- 
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ever, we can see that the averaged torque due to the disc that 
lies interior to the planet has converged to a value which 
is very close to that obtained in a laminar disc model. The 
outer disc torque, however, remains somewhat lower in am- 
plitude, such that the net torque experienced by the planet is 
positive over the simulation run time. It is well known that 
the stochastic nature of the disc turbulence, and the associ- 
ated stresses it generates, can lead to structuring of the disc 
radial density profile due to local variations in the effective 
viscous stress. This may provide a means of modifying the 
balance between inner and outer disc torques such that con- 
vergence to the expected type 1 value does not happen in 
practice for an individual planet in a disc, although it should 
happen on average for an ensemble of planets. One poten- 
tial way in which the underlying type 1 torques may not be 
recovered is if local variations in the viscous stress lead to 
undulations in the radial density profile. In this situation the 
corotation torque may become important in regions where 
the disc surface density increases locally with radius. The 
effect of this is illustrated by equation (|6]l, which shows that 
if the surface density profile becomes positive then the coro- 
tation torque experiences a boost because a becomes neg- 
ative. This phenomenon has been dubbed a 'planet trap' by 
Over longer viscous time scales such undulations in the 
density distribution will be smoothed out, but it is possible 
that planets may be able to migrate from one such density 
feature to another such that the next migration rate is sub- 
stantially reduced. 

Significantly longer simulations are required in order to 
examine the long term evolution of low mass planets em- 
bedded in stratified, global, turbulent disc models and to de- 
termine the efficacy of disc turbulence in modifying type 1 
migration. Furthermore, in order to build up a statistically 
significant picture of the distribution of outcomes, a number 
of such simulations will be required. In addition, the issue of 
how dead zones affect the stochastic torques also needs to be 
examined. It is worth noting that the simulation presented in 
this section required approximately 55 days of wall-clock 
run time on 256 processing cores. Clearly a full exploration 
of this problem remains a very challenging computational 
task. 

5.2 High mass protoplanets in turbulent discs 

The interaction of high mass, gap forming planets with tur- 
bulent protoplanetary discs has been considered in a number 
of papers (l61tl62t dTU IOI). In all these works, the disc model 
considered was non stratified such that the vertical compo- 
nent of gravity was neglected. It is well-known, however, 
that the flow of material in the vicinity of a giant planet de- 
velops a well-defined three dimensional structure, such that 
full 3D simulations are required to examine the details of the 
disc -planet interaction. 




Fig. 15 This figure shows a snapshot of the disc midplane density for 
a protoplanet with initial mass of 1/3 Mjup embedded in a turbulent 
disc. This snapshot corresponds to a time when the planet mass has 
reached ~ 4 Mjup by accreting gas from the disc. 



Figure [15] shows a snapshot of the midplane density for 
a turbulent disc with an embedded 4 Mjup protoplanet. The 
disc model in this case is the same as that described in sec- 
tion |5T| with H/R — 0.1. The planet had an initial mass 
Mp = 1/3 Mjup and was able to accrete gas from the disc. 
The gas accretion algorithm consists of the standard proce- 
dure of removing a fraction of the gas contained within one 
quarter of the planet Hill sphere at each time step. When 
magnetic fields are present in the simulation this procedure 
can lead to the formation of a magnetically dominant enve- 
lope forming around the planet. This problem was circum- 
vented by allowing the gas within the planet Hill sphere to 
have a finite resistivity such that the field diffuses from this 
region as it builds up. 

An initial planet mass of Mp = 1/3 Mjup is below the gap 
opening mass according to the criterion expresed by equa- 
tion ( fTsT l. As the planet accretes mass, however, gas is re- 
moved from the vicinity of the planet, and the growth of 
the planet mass allows tides to begin to open a gap once 
the Hill sphere radius begins to exceed the local disc scale 
height. This is in agreement with , who considered the 
transition from fully embedded to gap forming planets us- 
ing local shearing box and global simulations of turbulent 
discs. These simulations showed that gap formation begins 
when the disc reponse to the planet gravity starts to become 
non linear The combined Maxwell and Ryenolds stresses 
in simulations with zero-net magnetic flux (as considered 
here) typically give rises to an effective viscous stress pa- 
rameter ttvisc — 5 X 10^^, such that the viscous criterion for 
gap formation is satisfied when the criterion for non linear 
disc response is satisfied. Here avisc is defined to be ttvisc = 
(r) / (P), where T is the total (Maxwell plus Reynolds) stress. 
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Fig. 16 The left panel shows the gas density at the disc midplane, and the right panel shows contours of the quantity at the disc midplane. 



P is the thermal pressure, and angled brackets represent ver- 
tical and azimuthal averages. 

Global simulations allow the net torque on the planet due 
to the disc to be calculated and hence the migration time to 
be estimated. Simulations presented in (61; 10) for massive 
planets indicate migration times of ^ 10^ yr, in line with 
expections for type II migration. The left panel of figure [T6l 
shows the torques per unit mass experienced by the planet 
due to the disc as a function of time. Although gap forma- 
tion has not reached a state of completion (this is found to 
require approximately 400 planet orbits in laminar disc runs 
such as those presented in section |4]i, such that the torques 
have not reached their final values, this figure can be used 
to estimate a lower limit on the migration time. We express 
this as Tmig — j/{dj/dt), where j is the specific angular mo- 
mentum of the planet. The specific torque at f = 100 orbits 
is dj/dt = 2 X 10^^, (expressed in code units). The planet 
specific angular momentum j = \/3 such that the migration 
time corresponds to Tmig = 3 x lO"* years. Completion of gap 
formation will reduce the torques experienced by the planet 
considerably as the gas is pushed away from its vicinity by 
tides, causing the final migration time to significantly ex- 
ceed this value. But, as noted in section ISTl the run times 
associated with these full 3D stratified disc models only al- 
low full 3D magnetised and stratified disc simulations to run 
for ~ 100 orbits at the present time. 

A number of interesting features arise in simulations of 
gap forming planets embedded in turbulent discs, which dif- 
ferentiate them from runs performed using laminar but vis- 
cous disc models. Comparing figurefTSlwith figure|9]shows 
that the spiral arms in the turbulent disc appear less well- 
defined and more diffused than in the laminar models. This 
point is reinforced by the left panel of figure [17] which con- 
tains a snapshot of the planet-disc system in the r-(j) plane. 
The magnetic field topology is significantly modified in the 
vicinity of the protoplanet, as illustrated by the right panel 
in figure [TT] which displays a snapshot of the quantity in 
the disc in the vicinity of the planet. It can be seen that the 



magnetic field is compressed and ordered in the postshock 
region associated with the spiral wakes, increasing the mag- 
netic stresses there. Accretion of gas into the protoplanet 
Hill sphere causes advection of field into the circumplan- 
etary disc that forms there, creating a strong magnetic field 
in the circumplanetary disc, and this field links between the 
protoplanetary disc and the circumplanetary disc, producing 
an effect of magnetic braking of the circumplanetary mate- 
rial. Indeed, a direct comparison between a simulation of a 
3 Mjup protoplanet in a viscous laminar disc and an equiva- 
lent turbulent disc simulation presented by (10) suggests that 
mass accretion onto the planet in a magnetised disc may be 
enhanced by this effect. Simulations are underway at present 
to explore this effect within the 3D stratified model pre- 
sented in this section. It is worth noting, however, that these 
global simulations are of modest resolution, especially when 
considering complex physical processes occuring within the 
planet Hill sphere, and more high resolution work needs to 
be done to examine these issues in greater detail. 



6 Summary 

We have reviewed recent progress in the field of disc-planet 
interactions in the context of orbital migration. This progress 
has come about to a large degree from large scale two and 
three dimensional simulations that have utilised the most up 
to date supercomputer resources. These have allowed very 
high resolution simulations to be performed which exam- 
ine the detailed behaviour of the disc in the corotation re- 
gion, leading to a fuller understanding of the potential in- 
fluence of the corotation torque in modifying the migration 
of low mass and fully embedded planets. Simulations have 
also enabled the study of migration in non isothermal discs 
in which radiation transport is included, showing the impor- 
tance of relaxing the widely used assumption that the disc 
gas can be treated as being locally isothermal. The inclu- 
sion of MHD effects in full 3D global simulations is now 



15 




Time (Orbits) Time (Orbits) 

Fig. 17 The left panel shows the torque experienced hy the giant as a function of time. The right panel shows the evolution of the giant planet 
mass as it accretes gas from the disc. 



possible, allowing the influence of magnetohydrodynamic 
turbulence on planet formation and migration to be studied. 

Key questions for the future relate to how these various 
physical effects combine to change our current understand- 
ing of planetary migration. For example, the long-term ac- 
tion of corotation torques requires viscous diffusion to op- 
erate in the disc in order to prevent saturation. An important 
question is how these corotation torques operate in discs 
where the effective viscous stresses are provided by mag- 
netohydrodynamic turbulence ? Can the corotation torques 
remain unsaturated within a dead zone ? These and other im- 
portant questions are the subject of on-going research. The 
increasing power of available computing facilities will con- 
tinue to facilitate the development of increasingly realistic 
disc models, which will in turn lead to an increase in our 
understanding of planetary formation and migration. 
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